Preperation database to extract and stack asc files
A note to anyone who might happen to stumble across this… I am a beginner in R and have had no exposure to similar languages. I don’t know what I’m doing. The code herein is unlikely to be elegant and there are probably more efficient ways of running the code.
Built with ‘r getRversion()’.
You can load them using the following code which uses a function called ipak. Note this function checks to see if the packages are installed first.
packages <- c("tidyr", "rJava", "dismo")
source("../src/ipak.R")
ipak(packages)tidyr rJava dismo
TRUE TRUE TRUE suppressMessages(ipak)function(pkg){
new.pkg <- pkg[!(pkg %in% installed.packages()[, "Package"])]
if (length(new.pkg))
install.packages(new.pkg, dependencies = TRUE)
sapply(pkg, require, character.only = TRUE)
}write.csv(asclayerall, "../data/env/asc_stacklist.csv", header = TRUE)
Error in write.table(asclayerall, "../data/env/asc_stacklist.csv", header = TRUE, :
unused argument (header = TRUE)#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 1998)
ascname <- list.files("../data/env/ascii/ascii1998/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii1998/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist1998.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCLAYER ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 1999)
ascname <- list.files("../data/env/ascii/ascii1999/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii1999/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist1999.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2000)
ascname <- list.files("../data/env/ascii/ascii2000/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2000/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2000.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2001)
ascname <- list.files("../data/env/ascii/ascii2001/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2001/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2001.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2002)
ascname <- list.files("../data/env/ascii/ascii2002/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2002/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2002.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2003)
ascname <- list.files("../data/env/ascii/ascii2003/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2003/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2003.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2004)
ascname <- list.files("../data/env/ascii/ascii2004/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2004/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2004.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2005)
ascname <- list.files("../data/env/ascii/ascii2005/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2005/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2005.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2006)
ascname <- list.files("../data/env/ascii/ascii2006/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2006/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2006.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2007)
ascname <- list.files("../data/env/ascii/ascii2007/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2007/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2007.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2008)
ascname <- list.files("../data/env/ascii/ascii2008/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2008/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2008.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2009)
ascname <- list.files("../data/env/ascii/ascii2009/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2009/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2009.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2010)
ascname <- list.files("../data/env/ascii/ascii2010/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2010/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2010.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2011)
ascname <- list.files("../data/env/ascii/ascii2011/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2011/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2011.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2012)
ascname <- list.files("../data/env/ascii/ascii2012/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2012/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2012.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2013)
ascname <- list.files("../data/env/ascii/ascii2013/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2013/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2013.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2014)
ascname <- list.files("../data/env/ascii/ascii2014/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2014/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2014.csv", row.names = TRUE)
}#!!!!!! CHANGE ASCNAME, ASCPATH AND CSV FILE NAME
asclayer <- subset(asclayerall, ascyears == 2015)
ascname <- list.files("../data/env/ascii/ascii2015/", pattern = "*.asc", full.names = FALSE) # this folder currently contains other dummy files. its just a test to automate stacks
ascpath <- list.files("../data/env/ascii/ascii2015/", pattern = "*.asc", full.names = TRUE) # this folder currently contains other dummy files. its just a test to automate stacks
tempvar <- "temp.nc"
chlvar <- "chl.nc"
o2var <- "o2.nc"
salvar <- "salinity.nc"
mlpvar <- "mlp.nc"
sshvar <- "ssh.nc"
naosamplevar <- "naosample"
amosamplevar <- "amosample"
naoprevvar <- "naoprev"
amoprevvar <- "amoprev"
naowintervar <- "naowinter"
amowintervar <- "amowinter"
surfacevar <- "0.50576"
for (i in 1:nrow(asclayer)){
print(asclayer$stackname[i])
year <- asclayer$ascyears[i]
month <- asclayer$ascmonths[i]
depth <- asclayer$ascdepths[i]
for (k in 1:length(ascname)){
ascyear <- (sapply(strsplit(ascname[k], "_"), "[[", 1))
ascmonth <- (sapply(strsplit(ascname[k], "_"), "[[", 2))
ascvar <- (sapply(strsplit(ascname[k], "_"), "[[", 3))
ascdepth <- (sapply(strsplit(ascname[k], "_"), "[[", 4))
ascdepth <- gsub(".asc", "", ascdepth)
ascfile <- ascpath[k]
if (ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == depth){
asclayer$temp_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == tempvar & ascdepth == surfacevar){
asclayer$temp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == depth){
asclayer$salinity_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == salvar & ascdepth == surfacevar){
asclayer$salinity_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == depth){
asclayer$chl_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == chlvar & ascdepth == surfacevar){
asclayer$chl_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == depth){
asclayer$o2_depth[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == o2var & ascdepth == surfacevar){
asclayer$o2_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == mlpvar){
asclayer$mlp_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == sshvar){
asclayer$ssh_surface[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naosamplevar){
asclayer$nao_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amosamplevar){
asclayer$amo_sample[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == naoprevvar){
asclayer$nao_prev[i] <- ascfile
} else if(ascyear == year & ascmonth == month & ascvar == amoprevvar){
asclayer$amo_prev[i] <- ascfile
} else if(ascyear == year & ascvar == naowintervar){
asclayer$nao_winter[i] <- ascfile
} else if(ascyear == year & ascvar == amowintervar){
asclayer$amo_winter[i] <- ascfile
}
}
write.csv(asclayer, "../data/env/asc_layerlist2015.csv", row.names = TRUE)
}it may be better to have ascii lists per month (since models are per month)
split(as, as$ascmonths)
$`1`
$`2`
$`3`
$`4`
$`5`
$`6`
$`7`
$`8`